arXiv: 1509.05795vl [physics.class-ph] 4 Sep 2015 


Gain and noise spectral density in an electronic parametric amplifier with 

added white noise 

Adriano A. Batista^ and A. A. Lisboa de Souza^ 

^ Departamento de Fisica, Universidade Federal de Campina Grande 
Campina Grande-PB, CEP: 58109-970, Brazil 
^ Departamento de Engenharia Eletrica, Universidade Federal da Paraiba 
Jodo Pessoa-PB, CEP: 58.051-970, Brazil 
(Dated: September 22, 2015) 

Abstract 

In this paper, we diseuss the behavior of a linear elassieal parametrie amplifier (PA) in the presenee 
of white noise and give theoretieal estimates of the noise speetral density based on approximate Green’s 
funetions obtained by using averaging teehniques. Furthermore, we give analytieal estimates for parametrie 
ampliheation bandwidth of the amplifier and for fhe noisy preeursors fo insfabilify. To validafe our fheory 
we eompare fhe analytieal resulfs wifh experimenfal dafa obfained in an analog eireuif. We deseribe fhe 
implemenfafion defails and fhe sefup used in fhe experimenfal sfudy of fhe amplifier. Near fhe fhreshold 
fo fhe firsf paramefrie insfabilify, and in degenerafe-mode amplifieafion, fhe PA aehieved very high gains 
in a very narrow bandwidfh eenfered on ifs resonanee frequeney. In quasi-degenerafe mode amplifiealion, 
we obfained lower values of gain, buf wifh a wider bandwidfh fhaf is funable. The experimenfal dafa were 
aeeurafely deseribed by fhe prediefions of fhe model. Moreover, we nofieed speefral eomponenfs in fhe 
oufpuf signal of fhe amplifier whieh are due fo noise preeursors of insfabilify. The posifion, widfh, and 
magnifude of fhese eomponenfs are in agreemenf wifh fhe noise speefral densify obfained by fhe fheory 
proposed here. 
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I. INTRODUCTION 


Parametrically-driven systems and parametric resonance occur in many different physical sys¬ 
tems, ranging from the mechanical domain to the electronic, microwave, electromechanic, op- 
tomechanic, and quantum domains. In the mechanical domain we have Faraday waves [[U, in¬ 
verted pendulum stabilization, stability of boats, balloons, and parachutes [[21. A comprehensive 
review of applications in electronics and microwave cavities spanning from the early twentieth 
century up to 1960 can be found in Ref. (Sl. A few relevant recent applications, in micro and nano 
systems, include quadrupole ion guides and ion traps DU, linear ion crystals in linear Paul traps 
designed as prototype systems for the implementation of quantum computing [l5]-17]|, magnetic 
resonance force microscopy DU, tapping-mode force microscopy DU, axially-loaded microelec¬ 
tromechanical systems (MEMS) DTOlfTTI . torsional MEMS [fT^ . In the quantum domain we could 
mention wideband superconducting parametric amplifiers lfT3ll . squeezing in optomechanical cav¬ 
ities below the zero-point motion liT4ll . and parametric amplification in Josephson junctions DTSll 
Parametric pumping has had many applications in the field of MEMS, which have been used pri¬ 
marily as accelerometers, for measuring small forces and as ultrasensitive mass detectors since the 
mid 80’s DUl. An enhancement to the detection techniques in MEMS was developed by Rugar 
and Griitter IfTTlI in the early 90’s that uses mechanical parametric amplification (before transduc¬ 
tion) to improve the sensitivity of measurements. This amplification method works by driving the 
parametrically-driven resonator on the verge of parametric unstable zones. 

Here, we study a classical parametric amplifier both in theory and in experiment with analog 
electronics. We investigate signal and idler responses near the onset of the first parametric instabil¬ 
ity. We find the experimental results of gain in the signal and idler responses accurately described 
by the theory. We also investigate the noisy precursors of instability. We provide an analytical 
expression for the noise spectral density of a parametrically-driven oscillator with added noise. 
Again, mostly, we have very good agreement between experiment and theoretical predictions. 
The novelty compared with the work of Wiesenfeld et al. in the 80’s [fT8l [T9l is that we provide 
a simpler theoretical derivation for the noise spectral density (NSD) based on Green’s functions 
and the averaging method, in which the contribution of both Eloquet multipliers involved are taken 
into account. Eurthermore, we provide an analytical expression for the noise spectral density in 
terms of the parameters of the system and not in terms of the largest Eloquet multiplier, which 
is left unspecified in their theory. The analytical calculation of the Eloquet multipliers can still 


2 


be a daunting problem to overeome. Another differenee, is that below threshold there is no pe- 
riodie solution, the system is in quieseent mode, a system that was not treated by Wiesenfeld et 
ah. Beeause we add noise to a linear parametrically-driven system below threshold, we can fit the 
NSD when there is no pumping with the NSD of a harmonic oscillator with added noise and thus 
calibrate the noise level. If we had used Wiesenfeld’s theory the measured noise level would be 
off by 6dB. In their model there is always a parametric pumping, since the noise is treated as a 
perturbation around a deterministic limit cycle of a nonlinear dynamical system. Hence, the pos¬ 
sibility of calibration of the NSD around a harmonic oscillator limit is not possible, or at least is 
not clear, in their theory. 

The contents of this paper are organized as follows. In Sec. Q we present our theoretical 
model, in Sec. ( [1111 ) we describe our experimental setup, in Sec. ( [TV] ) we present and discuss our 
numerical and experimental results, and in Sec. (|V]) we draw our conclusions. 


II. THEORY 


The block diagram shown in Fig. is an schematic implementation of an electronic circuit of 
the parametric amplifier. The block with the x symbol in it represents a multiplier and the box 
with a S symbol is a summer. The equation that describes this system is given by 


ds^ ds 


AVp cos{4:7i f s)Vc + Vq cos(27r/sS -f 0o), 


( 1 ) 


where Vc is the voltage on the capacitor, L is the inductance, C is the capacitance, R is the re¬ 
sistance, Vp is the pump amplitude, A is the amplification factor of the multiplier (with units of 
volt~^), Vo is the signal voltage amplitude, 2/ is the pump frequency, fs is the frequency of the 
signal, and 0o is the phase of the signal. We can simplify Eq. ([T]) with the adimensional time 
t = uqs, where uq = 1/VLC, 7 = 1/Q = R/{uoL), Fp = AVp, Fq = Vq. We find 

dj^x dx 

— + x = -1-^ + Fp cos{2ut)x -f Fq cos{uA + 0o), (2) 

where x{t) = Vc{t), u = 27r//a;o, and Ug = 2Trfs/uo. 
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A. The Green’s function of the parametric oscillator 


The equation for the Green’s function of the parametrically-driven oscillator described in Eq. Q 
obeys 

^ + l + 7^-i^pCos(2a;f) 

We use the retarded Green’s function which is = 0 for t < t'. By integrating the above 

equation near t = t', we obtain the initial conditions when t = t' + O'*", which are G(f, t') = 0 
and = 1.0. Assuming the parameters = 0{e), with 0 < e << 1, we can write 

the Green’s function in the stable zone of the parametrically-driven oscillator in the first-order 
averaging approximation lilOl as 


G{t,t') = 6{t-t'). 


(3) 


with 


GnM = 


G{t, t-r) Ri Go(r) -f Gp{t, r), 

e-'rr/2 


(4) 


UJ 


coshfK r) sinfcer) -I — sinhi^ r) cosfcur) 


Gpit.r) = ——e sinh(K r) cos(a;(2f — r)). 

CJK 

where T = t — t'>0, k = — 5^, (3 = —Fp/Au, and 6 = VL/2u:, and fl = 1 




B. The ac-signal response in the parametric oscillator 


Here we present the theory on classical linear parametric amplification near the onset of the first 
instability zone of the parametric amplifier, i.e. we analyze the response of a parametric oscillator 
to an added input ac signal. In the following, we will use the Green’s function of Eq. Q to obtain 
the solution of Eq. This is given by 


x{t) = / G{t,t')Focos{uJst' + (j)o) dt', 


(5) 


since we assume the pump and the input signal have been turned on for a long time, any homoge¬ 
neous solution has already decayed. The signal response is given by 

[ dt' Go{t — t') cosiojst' + 0o) = Re / dt' Go{t — 


= Re 

= Re 






' —oo 
roo 


Gn(T)e 


lUJsT 


= Re 




( 6 ) 
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In order to obtain 


^0(0;,) = / 


1 

^ Jo 


g-7i/2g*‘^si 


cosh(Kt) sin(a;t) H— sinh(Kt) cos(a;t) 

K 


dt, 


the following integrals were utilized 


(7) 


) 

e[-^/2+*(‘^^±‘^)]*sinh(Kt)dt = 


7 — 2i{u)s + u) 


2 [f /2 — K — i{u}s ± ce)] [7/2 + K — i{us ± ce)] ’ 

K, 

[7/2 — K — i{us ± a;)][7/2 + k, — i{us ± ce)] ’ 


and we obtained 


GoK) « i 

OJ 2l 


7 - 2i{ujs + oj) 


2 [ 7/2 - K - i{oJs + a;)][ 7/2 + k - i{ojs + ce)] 
7 — 2i{ujs — u) 

2[^/2 — K — ii^LJs — a;)][ 7/2 + k — i{oJs — 1 ^)] 

5 r 1 

2 Yl/2 — K — i{(jJs + u})]Yi/2 + K — iiujs + u})] 

1 

[ 7/2 — K — i{ojs — ce)] [ 7/2 + K — i{us — ce)] 


( 8 ) 


Henee, with the help of the appendix ealeulations, we find the stationary solution of Eq. @ to 
be approximately 


x{t) 


|Re 


1 

h* 

n 

-1 

1 


2co; L 


gi[(2a;+aJs)tH-(/)o] 


2ce . [7/2 — ^ + ^s)] [7/2 + K + iix ) + tUs)] 


^i[{2u)-ujs)t-4ifi\ 


[7/2 — K + i{u — CJs)] [7/2 + K + i{u — Co’s)] 


(9) 


With Us ^ u and when the PA is pumped near the onset of the first instability, we find the idler 
response to be 


f dt' Gp{t, t') cos{ust' + 0o) ~ 

' —00 


^i[{2uj-LLJs)t-(j>o] 


2u \ [7/2 — K + i{u — Co;s)][7/2 + K + i{u — Co’s)] 

We also find that the envelope of the time series x{t) given in Eq. Q is approximately 

^i[{uj-uJs)t-4>o] 




2Co’ [7/2 — K + i{u — Co’s)] [7/2 + K + i{u — Co’s)] 


( 10 ) 
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C. Noise spectral density in the parametric oscillator 


We will now investigate the effeet of noise on the parametrieally-driven oseillator [|2T]| . The 
parametric oscillator equation with added noise is given by 


X = —X — 72 ; + Fp cos{2ut) X + r{t), 


( 11 ) 


where r(t) is a Gaussian white noise that satisfies the statistical averages (r(t)) = 0 and 
(r(t)r(t')) = 2D5{t — t'), where D is the noise intensity. We shall now use the same analyt¬ 
ical method developed in Refs. [|2T1I^ to study the parametrically-driven oscillator with noise as 


given by Eq. ( [TT] ). 

Using the Green’s function we obtain the solution x{t) of Eq. ([IT]), which is given by 


x{t) = / dt' G{t,t')r{t'), 


( 12 ) 


assuming that the added noise has been turned on for a long time, hence, any homogeneous solu¬ 
tion has already faded out. 

The Eourier transform of x{t) is given by 


x[u = 


F''^x{t)dt = / dt'rit') / dt [Go(t — t') + Gp(t, t')] 


(13) 


= r{iy)Go{iy) + / dt'r{t') / dt F''^Gp{t,t') 


Hence, with the help of the appendix calculations, we find 

f(z/ -f 2uj) ^ f(z/ — 2uj) 

[p- + i{iy + cj)] [p+ + i{iy -f ce)] [p_ + i{iy - ce)] [p+ + i{u - u)]\ ’ 

(14) 

where p± = — 7/2 ± k are the Eloquet exponents [l20l in the first-order averaging approximation. 
Since the stochastic process of a parametrically-driven oscillator with added noise as defined in 
Eq. ( [TT] ) is cyclo-stationary, the correlation function of x{t) is not translationally invariant and the 
Wiener-Khinchin theorem is not valid. Therefore, one needs to perform a time average over the 
usual NSD [fT9ll . One then finds 


x(z/) = r{iy)Go{iy) - ^ 


/ CO roo 

dTe~^'^'^{x{tT)x{t)) = — du'e'^^’-'~'^'')^{x{—u)x{u')) 

2Trl_^ 


= lim 


Ai/-s>0+ 




27r 


(15) 
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With the help of the relation (r(z/)r(z/')) = 4:TtD6{i> + u'), we find, where k = ik!' is imaginary 

r 1 


S^{u)^2D{\G^{u)\^ + 


[72/4 -\- (^ i ! ijj — [^2^4 _j_ _|_ /-")2j 

1 


[72/ A+ {v — u — k"Y] [72 /4 + (i/ — ce + k")2] 

or where k is real (nearer the onset of instability), 


(16) 


S^{u)^2D{\Go{y)\^ + 

+ 


4ci;2 [[(7/2 - k)2 + (i/ + a;)2][(7/2 + k)2 + (z/ + u : Y ] 
1 


(17) 


(18) 


[(7/2 - k)2 + (z/ - Ci;)2] [(7/2 + k) 2 + (z/ - Ci;)2] _ 

When the pumping is turned off {Fp = 0) in Eq. ( [T4| ) and we take cu = 1, we obtain 

2[72/4+(z/-1)2]’ 

whieh is a very good approximation of the harmonie oseillator noise speetral density in high-Q 
oseillators. Near the instability threshold and with u Ri u, we obtain 

1 —7 + 2i[5 + V — u] 


^ dice [7/2 — K — i(z/— ce)] [7/2 + K — i(z/— cu)] ’ 

and the NSD is given approximately by 


(19) 


S^{v)^2D{\Go{i^)\^ + 
D 
D 


4a;2[(7/2 — k)2 + (z/ — a;)2] [(7/2 + kY + (z/ — a;)2] 
7^/4 + [5 + z/ - cu]^ + 


2a;2[(7/2 — kY + [v — a;)2] [(7/2 + kY + {v — uY ] 
7^/4 + (5 + z/- cj)^ +1 


Au'^jk (7/2 —/s;)2 + (z/— a;)2 ’ 

whieh is not exaetly a Lorenzian eurve as predieted by Wiesenfeld et al. 


( 20 ) 


III. APPARATUS 

In this seetion, we deseribe the eleetronie eireuit eoneeived to implement a parametrie amplifier, 
whieh is shown in Fig.|^ 

A. Analog electronic circuit of the parametric amplifier 

The eore of the parametrie amplifier is shown in Fig.|^and is eomprised of 3 main eomponents: 
a 4-quadrant analog multiplier (AD633), whieh has a eonversion gain of 1/10F“^, a (weighted 


7 

















inverting) summer implemented by one operational amplifier (opamp), and an eleetronie induetor 
(IH) implemented with two opamps (the well-known Antoniou induetor-simulation eireuit). Un¬ 
like the parametrie oseillator eireuit of Ref. [|2^ . we use a linear eapaeitor in plaee of the nonlinear 
varieap diode, whieh makes our eireuit simpler to analyze. Although nonlinear behavior will even¬ 
tually appear onee the threshold to instability is erossed, we are interested in the operation of the 
eireuit as an amplifier in the linear regime, and not as a nonlinear oseillator. 

With the ehoiee of C 2 = 552pF, a nominal resonanee frequeney of 677AHz is obtained. When 
ehoosing the eapaeitor C 2 , eare must be taken to avoid lowering the quality faetor of the resonator 
(whieh is otherwise determined by the quality faetor of the induetor along with the value of the 
resistor i?10). 

The dynamieal variable of the eireuit is the voltage at the eapaeitor C 2 , whose readout is 
buffered before being eonneeted to instruments to avoid eurrent drains during measurements. To 
avoid outside eleetromagnetie interferenee, the PA eireuit was enelosed in a metallie box. 


B. Measurement Setup 

The setup used to eharaeterize the behavior of the parametrie amplifier is shown in Fig. 

It is eomposed of a veetor signal analyzer (model HP89A10A), a waveform generator (model 
33210A), and an oseilloseope (model MSO — X202AA). The waveform generator is used to 
generate the (sinusoidal) pump signal. The veetor signal analyzer has an internal souree that is 
used as the signal to be amplified. Moreover, the analyzer has two input ehannels, whieh ean be 
eonfigured to evaluate, for example, speetrum ratios. Besides the frequeney-domain analysis with 
the HP89A10A, we have also observed signals in time-domain with the aid of the oseilloseope. 

The setup is flexible enough to enable the experimental eharaeterization of behavior of the Para¬ 
metrie Amplifier under different operating eonditions, depending upon whether a short-eireuit or a 
signal souree is eonneeted to the inputs of the amplifier. The different measurement eonfigurations 
are deseribed below. 
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IV. NUMERICAL AND EXPERIMENTAL RESULTS 


A. Harmonic oscillator resonant curve 

To evaluate the response of the harmonic oscillator, for which Fp = 0 in Eq. Q, the pump input 
port of the amplifier is short-circuited, while the signal input is fed with a signal source (please 
refer to Fig. connections b and d). We have measured the response (gain) of the harmonic 
oscillator for a broad range of frequencies. The results are shown in Fig. One can see that the 
equivalent quality factor Q of the circuit was found to be about 65. This means that if we consider 
the circuit of Fig. the equivalent series resistance is about 6190. In addition to the 560 shown in 
Fig.[^ there are additional losses attributed to the electronic inductor, since the quality factor of the 
capacitor C 2 was independently measured as about 124. The resonance frequency (/o = uiq/ (2vr)) 
was found to be around 6400if^, a value less than 6% lower than the nominal frequency. 


B. Instability boundary 

After obtaining the resonance curve of the circuit in harmonic oscillator configuration, we have 
evaluated the instability boundary of the circuit in parametric oscillator configuration. This con¬ 
figuration is obtained by pumping the amplifier while short-circuiting its signal input. In Fig. 
this corresponds to connections a and c. For each value of pump frequency set, the pump am¬ 
plitude was slowly increased from a very low initial value until an oscillation of more than 0.5F 
in amplitude was observed at the output of the amplifier. When comparing the experimental data 
against the numerical data from Eq. 0. we have found that the equivalent Q of the resonator is 
about 65 (see Fig. [^, a value consistent with the one observed for the harmonic oscillator resonant 
curve. Moreover, the instability line is centered on the same value of frequency found for the peak 
of harmonic oscillator resonant curve (about bdOOFT^), which represents half the value of pump 
frequency for which the lowest pump amplitude is necessary for instability. 

Hereafter the parameter values Q = 7“1 = 65 and ccq ~ 40212 rad/s will be used in fitting the 
experimental data of parametric amplification against theoretical predictions. 
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C. Parametric amplification 


Finally, we set the eireuit in parametrie amplifieation eonfiguration, by inputing both pump and 
signal, as deseribed in Fig.j^with eonneetions a and d. In all the results presented below, the input 
signal amplitude was limited to 2mV. Owing to the high gains whieh ean be obtained, this value 
has to be kept small enough to avoid saturation of the amplifier due to limitation of supply voltage 
bias and intrinsie nonlinearities of the aetive eomponents of the eireuit. 

In Fig. we show a time series obtained from the eireuit set with pump amplitude Vp = 31/, 
pump frequeney I. 8/0 and input signal frequeney fs = 0.95/o, along with the envelope predieted 
by the Eq. ( fTO] ) with Fp = 0.3. 

In Fig.[^ we show a time series obtained from the eireuit set with pump amplitude Vp = 0.29V, 
pump frequeney 2/o and input signal frequeney fs = 0.99/o, along with the envelope predieted 
by the Eq. ( fTO] ) with Fp = 0.029. The best agreement between the experimental time series and 
the theoretieal envelope is obtained under this eondition (degenerate-mode amplifieation). This 
oeeurs beeause the aeeuraey of the perturbative methods used (averaging and harmonie balanee) 
is higher the smaller the parameter Fp is. 

In Eig. we show a time series obtained from the eireuit set with pump amplitude Vp = 41/, 
pump frequeney 2.2 fo and input signal frequeney fs = 1.05/o, along with the envelope predieted 
by the Eq. ( [T0| ) with Fp = 0.4. 

In Eig. 1^ we show the signal speetrum measured at the output of the amplifier with the aid 
of the signal analyzer (please refer to Eig. |^. The experimental eonditions eorrespond respee- 
tively to those assoeiated with the time series presented in Eigs.|^|^ and[^ The signal analyzer 
eomputes the speetrum eoneurrently to the aequisition of the waveforms by the Oseilloseope. The 
experimental results are eompared against the numerieal Eourier transforms of the time series x{t) 
obtained from numerieal integration of Eq. 0. The position and magnitude of the main peaks 
(signal and idler) is in agreement with the theoretieal predietions of Eq. In frame (b), the 
middle peak in the experimental data is a noisy preeursor of the parametrie instability, whereas 
in frames (a) and (e) one hardly notiees any effeets of noise. To help visualize these effeets, the 
insets in frames (a) and (e) show details of noisy preeursors in the experimental data. Eurther be¬ 
low we will explain why these preeursors are more relevant in the degenerate-mode amplifieation 
of frame (b), than in eases (a) and (e). In all eases, though, the response of the PA to noise is 
more pronouneed when the amplifier is tuned to the vieinity of the transition line [l20l - l^ . Henee, 
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the closer one gets to the transition line, the higher the noisy precursor lines in the spectrum will 
be. A detailed analysis of the noise spectrum is made in the next subsection. The noisy precursor 
effect is in qualitative agreement with the observations of Jeffries and Wiesenfeld on the effect of 
broadband noise on the power spectrum of coherent signals, which was first investigated (theoreti¬ 
cally and experimentally) near period-doubling and Hopf bifurcations in a periodically driven p-n 
junction [fT^[T9l and in the context of parametric amplification in Josephson junctions Il24l . 


In Fig. 10 one can see the signal and idler responses of the PA as a function of frequency. 
Again, one can see that there is a much broader bandwidth of gain in quasi-degenerate mode of 
amplification, such as in frames (a) and (c), than in degenerate-mode amplification, frame (b). 
Moreover, the bandwidth of gain, the peak positions, and the gain xbandwidth product can be 
tuned as well, unlike in the degenerate-mode amplification. 


D. Noise spectral density 


The spectra of the signal responses can be used to predict where the lines due to noise will turn 
up in the Fourier transform of the time series x{t), as shown in Fig. Since there is always a 
small amount of added noise of very broad bandwidth, which comes along with the input ac signal 


or is intrinsic to the circuit, there will be noise components everywhere in the spectra of Fig. 10 


Hence, as our system is linear, the spectral components of noise will be amplified in the same 
way as the input ac signals are amplified. The strongest case is seen at the peak of degenerate¬ 
mode parametric amplification, when ojg = 1, where the noise will be amplified roughly by 55dB. 
The corresponding noise line can be seen in Fig. |^b). The noise lines in quasi-degenerate-mode 
parametric amplification are much smaller, since the the peaks in frames [T0]^a) and (c) have gains 
of only roughly 39dB and 42dB, respectively. Nonetheless, if one looks in the insets of frames 
(a) and (c) of Fig. one can see elevations in the noise level exactly at the peaks of the signal 
response of Fig. [T^ The horizontal and vertical dashed lines at the signal peaks of Fig. [T0|are 
reproduced in the insets of Fig. for clarity. One can then compare the magnitudes of the noise 
line peaks of Fig. The difference in gain at noise line peaks between frame |^b) and (a) is 
20 logio( 4 . 5 / 0 . 55 ) = 18.2dB, whereas between frame (b) and (c) is 20 logio( 4 . 5/1.1) = 12.2dB. 
On the other hand, the predicted differences in noise peaks based on theory for the signal gain 


shown in Fig. 10 is 16dB and 13dB, respectively. These results show that our experimental results 
are consistent with the predictions of our proposed linear theory to within an error of under 2.5dB. 
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We note that the same noise signatures appear in the speetrum when there is no input ae signal, 
eonfirming that the noise precursors are a consequence of linear response and a consequence of 
the signal gain spectra of the PA, such as the ones presented in Fig. [T^ 

In Fig. we show the noise spectral density for our circuit setup in harmonic oscillator mode. 
Here we fit the data with a noise intensity of D = 3.08 x 10“®V^/Hz. The source of this noise 
is intrinsic to the circuit. Here we showed the exact theoretical result for the NSD of a harmonic 
oscillator process given by 5'° alongside the approximate result Sx from Eq. ( [T8] ). In our units 
Sx has dimensions of since the time is adimensional. In order to set it in proper units, Sx is 
divided by uq. Both theoretical predictions yield nearly the same spectrum and account well for 
the experimental data. On the other hand, the NSD from Eq. (13) of Ref. [|T8]| gives a result that 
is 6 dB higher than the measured noise. Apart from this, their result can be rescaled to be exactly 


our result from Eq. ( [18] ) if one divides their prediction by 4. Here we have a way of measuring 
the noise level when there is no pumping, whereas in their model, there is no quiescent solution 
without noise. They developed a perturbative solution due to noise around a limit cycle (an isolated 
periodic solution) and did not calibrate their solution in the zero pumping limit. 

In Eig. [^one can see the comparison between theoretical predictions, given by Eq. ( [T7| ) or 
by Jeffries-Wiesenfeld model, and experimental results for the NSD of a parametrically-driven 
oscillator with added white noise. Here, the noise level used to fit the data was the same one of 
the harmonic oscillator configuration. In frames (a-c) and (e-g) the Eloquet exponents are still 
complex, as can be seen in the gap between the noise peaks. These peaks are located symmetri¬ 
cally with respect to oo and the distance between them is twice the imaginary part of the Eloquet 
exponents. In frame (d) the Eloquet exponents are real and, hence, there is no split peak in the 
noise spectrum. The discrepancy of theory and experiment at the peak of the NSD is certainly 
due to nonlinear effects not accounted for by our linear model. Again here the Jeffries-Wiesenfeld 
model is off with a level of noise higher by slightly over 6 dB. Note also that the absence of sharp 
peaks in the spectrum is because we have a quiescent solution and not a periodic orbit solution 
when there is no noise. 


V, CONCLUSION 

Here we obtained experimental results on gain and bandwidth of classical parametric ampli¬ 
fication that are quantitatively well approximated by a theory based on averaging techniques and 
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on Green’s funetion theory. Although one ean reaeh extremely high values of gain in degenerate 
mode parametrie amplifieation, the eorresponding bandwidth is very narrow. This is sometimes 
an undesirable eharaeteristies, sinee it makes tuning to the signal a diffieult task. We have found 
that the bandwidth of the parametrie amplifier ean be inereased if we set the amplifier in quasi- 
degenerate-mode amplifieation, i.e. with detuning (fl 7 ^ 0). This eomes at the expense of the high 
gain obtained in degenerate-mode amplifieation. Guided by the model developed here, the opti¬ 
mal amount of gain and bandwidth eould be found by earefully tuning the pump frequeney and the 
pump amplitude. Moreover, the bandwidth of gain, the peak positions, and the gain xbandwidth 
produet ean be tuned as well, unlike in the degenerate-mode amplifieation. 

Furthermore, we presented a theory for obtaining an approximate analytieal expression for the 
NSD of the parametrie oseillator with added noise that aeeounts for the noisy preeursors of insta¬ 
bility in the PA. With the information of the signal and idler gain speetra and the input noise level, 
one eould, in prineiple, determine how the noise features will be manifested in the power speetrum 
of the output signal of the PA. Although Wiesenfeld et al. [fT^ IT9]| developed a general eompre- 
hensive model for aeeounting for the noisy preeursors of bifureations of eodimension -1 based on 
Floquet theory, the model we propose is simpler to apply and has all the theoretieal expressions 
neeessary to eompare with experimental results. It does not depend on generie unspeeified Floquet 
exponents. It is worthy to mention that the present model deseribes the simplest nontrivial system 
to whieh the generie theory of Wiesenfeld eould be applied to, but whieh was aetually overlooked 
in the literature so far to the authors knowledge. Furthermore, we have results for the noise speetral 
density that take into aeeount the presenee of two Floquet multipliers, and not just one as deseribed 
in Ref. [|T9l . The existenee of two eomplex Floquet multipliers is eharaeterized by the appearanee 
of two noise peaks in the NSD speetrum, whieh ean be seen in several of our results. This is 
espeeially relevant in high Q parametrie oseillators, where the stable parameter region in whieh 
the Floquet amplifiers are real deereases with inereasing Q. Henee, it beeomes harder to tune into 
this region, speeially, when there is detuning in the pump frequeney with respeet to the resonanee 
frequeney, so one has to aeeount for the eontribution from both Floquet exponents. Also, due to 
this, the noisy preeursors speetral eurves are not properly Lorenzian anymore. 

The predieted gain eurves proposed here ean be used to determine signal and idler gains, the 
output noise speetral density, and the figure of noise of the PA (ratio of output and input signal- 
to-noise ratios of the PA). Furthermore, our PA eireuit ean be used as a simple and inexpensive 
experimental platform to test reeent theoretieal predietions [l 20 l - l^ on thermal noise squeezing. 
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Finally, the close proximity of the theoretical predictions and the experimental results obtained 
here indicates that one could design electronic devices based on PAs that could achieve extremely 
high gains, have very little noise, and be tunable. Future experimental and theoretical work on the 
PAs will be performed in analyzing the effects of nonlinearity on the dynamic range of amplifica¬ 
tion and the phenomenon of noise squeezing. 


Appendix: Idler response calculations 


The idler response is given by 


dt' Gp{t, t') cos{uJst' + 0o) = / dr Gp{t, t — r) cos{uJs{t — r) - 1 - 0o) 


= —— I dr e smh{KT) cos{u{2t — t)) cos{us{t — t) + (l)o) 

dr e“''"^^^sinh(fi;r) [cos ((a; + Us){t — t) + ut + 0 o) + cos ((a; — Us){t — r) + ut — 0 o)] 


UK JO 

— / 

2uk Jo 

■—Re [ dr sinh(Kr) ^e^lG+‘^s)d-r)+<t>o] + gi[(a;-a;,)(t-r)- 0 o]j 

2uk I Jo 


= -ARe 




+ 


2u [^/2 — K i{u -f cj 5)][^/2 -|- k i{u -f t.i.’s)] 

^i[{2u)-uJs)t-4io\ 


[ 7/2 - K + i{u - Us)\Yl/2 + K + i{u - cjs)] J 

where we used the following result 

POO 1 

/ dr e-^G‘i+iG±^s)]T r) = - 


(A.l) 


7/2 — K + f(a; ± Ws) ^/2 + K + i{u ±Us) 

K 

[ 7/2 — K + i{u ± Us)]Yi/2 + K + i{u ± Us) 


(A.2) 


Appendix: Noise response 


We have to solve 



(A.l) 
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dt eI-T'/2+*(‘^+^)li sinh(K (t - t')) = h + l 


2 , 


where 


Ii = 


„-Kt' poo \-^/2+i(u+uj)\t' 

^ _ / Af A-^/2+i(u+u,)]t Kt ^ _ 

2 ^-2k-2i{i2 + u) 


^Kt' roo 


J 2 =-/ dt e^-'r/2+i(i'+u^)]t^-Kt ^ 


g[-7/2+i(i/+a;)]t' 

7 + 2/t; — 2i(z/ + uj) 


where 


Ii + I 2 — e 


— ^[-7/2+7^^+!^)]*' 


[7 — 2fi; — 2 i{uj + u) 'y + 2 k — 2 i{u + u) 
4,^e[-7/2+*(‘^+i")]t' 


K, = 


[y — 2 k — 2 i{u + z/)] [7 + 2 /t; — 2 i{u + z/)] 

poo 

/ dt e[-^/2+*('^-^)li sinh(K (t - t')) = K^ + K 2 , 

Jt' 

„-Kt' poo [_^/2+i(j^—a;)]t' 

— / dt el-i/2+.(--»)lie-.i = ^^^ 
2 oj — 2 k — 2i(iy — uj) 


-^Kt' /‘OO 


K2 = -/ dt e[-T'/2+*(^-‘^)]ie-'^‘ = 


g[- 7 / 2 +i(i/-a;)]t' 

■y + 2 k — 2 i{v — Lo) 


/ dte^‘'^G„(t,t') ^ -—e 

li-t KU 


+ 




{h + l2)e^'' + {Ki + ir2)e-'‘"''] /2 

^i(v’+2u)t' 

2(0 [[7/2 — K —-z(z/+ a;)][7/2 + K — ^(z/+ cj)] 
[7/2 — K — i(z/ — a;)][7/2 + k — liyu — a;)] 


dt'r(t') f (it e*^*Gp(t, t') ~ 

lit 


+ 


r(z/ + 2 (o) 


[7/2 - K-i{u + a;)][ 7/2 + k - i{u + a;)] 
f(z/ — 2 u) 


['y /2 — K — i{v — cj)] ['y /2 + K — i{v — a;)] 
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FIG. 1. Schematic block diagram of the proposed analog electronic implementation of the parametric 
amplifier. The box on the left represents a multiplier with gain A and the box on the right is a non-inverting 


summer. 



FIG. 2. Electronic circuit implementing the parametric amplifier. The IH electronic inductor is im¬ 
plemented with the aid of operational amplifiers. The IC AD633 represents the multiplier,with a gain of 
10E“^. For simplicity, the output buffer and the power supplies are not shown. 
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FIG. 3. Setup used to characterize the behavior of the parametric amplifier, allowing different operating 
conditions. With options h and d, the PA behaves as a driven harmonic oscillator. With options a and c, the 
PA behaves as a parametric oscillator. With options a and d, the PA behaves as a parametric amplifier if the 
pump is kept below the instability threshold. 
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fs (Hz) 


FIG. 4. Harmonic oscillator (pump off) resonant curve, obtained using connections b and d in Fig.[^ Here 
we fit the experimental resonant curve with the theoretic prediction and find a Q of approximafely 65 and a 
resonanf frequency of /o = 6400Hz. 
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Boundary of the 1 st instability zone 



FIG. 5. The instability boundary of the amplifier circuit was obtained with the circuit set in parametric 
oscillator mode, as illustrated in Fig.|^ Above the boundary line, the quiescent state of the circuit (I 4 = 0) 
is unstable and the circuit becomes a nonlinear oscillator, whereas below this line the state I 4 = 0 is stable. 
Very good agreement between experiment and theory was found with the quality factor Q = 65. 
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Q=65,(0=0.9, O)g=0.95, Fp=0.3 



t (sec) 


FIG. 6. Comparison between a time series from experimental data of 14(f) and corresponding envelope 
obtained from theory in Eq. ( [T0| ). The amplification here is set in quasi-degenerate mode with all parameters 
indicated in the figure. 
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FIG. 7. Comparison between time series of experimental data of 14(f) and envelope obtained from theory 
in Eq. The amplification here is set in degenerate mode with a small detuning between half the pump 
frequency, cj = 1, and the input signal frequency, Ug = 0.99. 
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Q=65,(0=l.l, (Og=1.05, Fp=0.4 



t (sec) 


FIG. 8. Comparison between a time series from experimental data of Vc{t) and envelope obtained from 
theory in Eq. ( [T0| ). The amplification here is set in quasi-degenerate mode. All parameters are indicated in 
the figure. 
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FIG. 9. Comparison between experimental data of Fourier transform of V"c(f)/Vo and numerical Fourier 
transform of the time series x{t) /Fq from Eq. in The Fourier transforms refer to the time series presented 
in previous figures. The position and magnitude of the main peaks (signal and idler) is in agreement with 
the theoretical prediction of Eq. The insets in frames (a) and (c) zoom in the noise signatures of the 
experimental Eourier transforms. In frame (b) the middle peak in the experimental data is a noisy precursor 
of the parametric instability. 
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FIG. 10. Theoretical signal and idler response gains compared with RLC resonant curve gain. From Eq. Q 

the signal gain in dB scale is given by 20 logig(|G'o(a;s)|), with Go{uJs) given in Eq. ([^. Similarly, the idler 

gain is obtained from Eq. Q. In frames a) and c) the Eloquet exponent is complex, what causes the split 
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peaks in the signal and idler responses, whereas in frame b) it is real. As a comparison, we show the gain 
of the harmonic oscillator response (dot-dashed line). 
























FIG. 11. Comparison between theoretical predictions and experimental results for the noise spectral density 
of a harmonic oscillator. Here Fq = 0. The measuring apparatus depicted in Fig. was setup such that the 
Signal and the Pump ports were short-circuited (options b and c). The results for Sx are given by Eq. ( [TT] ) 
without parametric pumping and at resonance, whereas is the exact harmonic oscillator spectral density. 
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Noise Spectral Density in dB 



FIG. 12. Comparison between theoretical predictions, given by Eq. ( [T7| ) and experimental results for the 
noise spectral density of a parametrically-driven oscillator with added white noise. Here Fq = 0. This was 
achieved in the setup by short circuiting the input signal (option c in Fig. [^. The noise level used here is 
the same one used in the harmonic oscillator driven by white noise of Fig.fTT] 
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